A Generalized Ginzburg-Landau Approach to Second Harmonic Generation 
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We develop a generalized Ginzburg-Landau theory for second harmonic generation (SHG) in 
magnets by expanding the free energy in terms of the order parameter in the magnetic phase and 
the susceptibility tensor in the corresponding high-temperature phase. The non-zero components of 
the SHG susceptibility in the ordered phase are derived from the symmetries of the susceptibility 
tensor in the high-temperature phase and the symmetry of the order parameter. In this derivation, 
the dependence of the SHG susceptibility on the order parameter follows naturally, and therefore its 
nonreciprocal optical properties. We examine this phenomenology for the magnetoelectric compound 
Cr203 as well as for the ferroelectromagnet YMnOs. 

PACS Numbers: 42.65.-k: Nonlinear Optics, 75.50.Ee: Antiferromagnetics 



Second Harmonic Generation (SHG) is a very useful 
technique to study the nonlinear optical properties ||l| in 
magnetic materials. The recent observation [plH] of non- 
reciprocal optical effects (i.e., not invariant under time re- 
versal operation) in the magnetoelectric material Cr203 
below the Neel temperature T/v has been of great im- 
portance in the study of antiferromagnetic (AFM) or- 
dering by light. With the help of SHG, photographs of 
the antiferromagnetic domains 0,0] in Cr203 were ob- 
tained what confirms that these experiments can distin- 
guish between the two magnetic states that are related to 
each other by the time reversal operation, and therefore 
its nonreciprocity. These observations can be explained 
by an interference effect between a time-symmetric mag- 
netic dipole contribution and a time-antisymmetric elec- 
tric dipole contribution 0. Soon after these experi- 
ments were done, a microscopic theory was proposed |5|] 
which could explain quantitatively the non-reciprocal ef- 
fects observed in Cr203. It was shown that the electric 
dipole contributions were linearly proportional to the an- 
tiferromagnetic order parameter giving rise to the time- 
antisymmetric character to the electric dipole tensor. 
The study of this dependence is of crucial importance in 
order to understand the non-reciprocal character of the 
SHG tensors. The derivation of a micoscopic theory for 
a specific effect, in our case SHG, can be rather compli- 
cated depending on the properties of the material under 
study ||g,g] . Concrete information of the SHG process in 
materials where a transition takes place, i.e. paramagnet- 
antiferromagnet or paraelectric- ferroelectric, can be ob- 
tained at a simpler level, i.e. by considering only the 
symmetry arguments. In particular, we are interested in 
investigating the dependence of the SHG susceptibilities 
in the ordered phase on the order parameter. 

A powerful phenomenological theory suitable to de- 
scribe phase transitions is the Ginzburg-Landau theory. 
The Ginzburg-Landau approach is based on the existence 
of an order parameter in the ordered phase and on sym- 
metry considerations m. Pershan in 1963 M showed 



that the tensors for nonlinear electro- and magneto-optic 
effects could be derived from a phenomenological time- 
averaged free energy. It is our purpose in this paper to 
extend the formulation of Pershan by including the order 
parameter explicitly in the expansion of the free energy 
for SHG. We shall show that the non-zero components 
of the SHG susceptibility tensor in the ordered phase are 
naturally obtained from the symmetry of the susceptibil- 
ity tensor in the high-temperature phase and the sym- 
metry of the order parameter. Once the dependence of 
the SHG susceptibility tensor on the order parameter is 
known, the non-reciprocal optical properties below the 
transition temperature follow. We explicitly verify this 
formulation in the magnetoelectric compound Cr203 as 
well as the hexagonal ferroelectric-antiferromagnetic ma- 
terial YMn03. 

A convenient starting point in order to describe the op- 
tical nonlinearities is given by a time-averaged free energy 
F as proposed by Pershan ^ . A dipole expansion of the 
induced current J = dP/dt + cV x M-9(V-Q)/at + . . ., 
where P, M and Q are respectively the electric dipole po- 
larization, the magnetization and the electric quadrupole 
polarization, leads to Q 

n 

F = -2i?e^ [E*(w,) • P{uj,) + H*{uj,) ■ M(w,)] , 

(1) 

where we have discarded the electric quadrupole and 
higher-order terms and a constant contribution. Here 
E(wi,), P{uj,y), il{uji,) and M(a;i,) are respectively the 
Fourier components of the electric field, electric dipole 
polarization, magnetic field and the magnetization and v 
denotes the number of partial waves. The usual relations 



P(c^.) 



-aF/9E*(cj^) 
-dF/dIl*{uj^) 



are fuUfiUed. The electric dipole polarization and the 
magnetization are, in general, nonlinear functions of E 



(and H) which may be expanded into a power series in 
E (and H) . The second order term for the electric dipole 
polarization is 

P,{2lo) = x\f^\2u,Lo,Lo) E,{lo) E^ico) (2) 

which corresponds to the SHG electric dipole contribu- 
tion (ED). The corresponding free energy is then (sec Eq. 
(i): 



p{ED) 



{ED) 



X^;,>{2^,^,u) E*{2lo) E,i^) EkH 
+X*J''''\2uj,uj,uj) E,{2lu) E*{oj) E^co)] (3) 



Summation over repeated indices is implicit in the above 
formulas. 

As a second example, we consider the magnetic dipole 
contribution (MD) to SHG, 

M,{2uj) = xl^r^Cac^.c^.c.) E,{cu) Ekico) , (4) 

leading to 



F 



{MD) _ 



x'^^l''\2uj,uj,uj) H*{2uj) E,{oj) Ek{oj) 



+x:J^'''\2cj,cj,cj)H.i2cj)E;icj)El{cj)] . (5) 



piMD) jg |.]_^g j?j.gg energy corresponding to the MD con- 
tribution to SHG. Both the ED and the MD contribute 
to the source term S, i.e., as ^o[(5^-P/3i^) + V x (dM/dt)] 
in the wave equation for the electric field, derived from 
the Maxwell's equations in classical electrodynamics Q. 
The measured output intensity in a nonlinear optical ex- 
periment is / ex I 5* I and once the SHG susceptibility 
in the ordered phase is known, possible non-reciprocal 
properties of the system may be derived. 

In the context of the antiferromagnetic Cr203 it has 
been observed ||] that xljk exists only in the ordered 
phase (T < Tjv) and it has been shown MM that it is lin- 
early proportional to the antiferromagnetic order param- 
eter, as it should be if we had defined a Ginzburg-Landau 
free energy for the ordered phase. In the framework of the 
classical Ginzburg-Landau approach to phase transitions 
the allowed terms contributing to the free energy follow 
from the symmetry of the order-parameter and from the 
symmetry of the lattice @,^. It is therefore of inter- 
est to ask oneself whether one can combine the standard 
Ginzburg-Landau approach for the magnetic properties 
of e.g. Gr203 with the expressions Eq. (|^) and Eq. (g) in 
order to obtain a more general formulation for the non- 
linear magneto-optical properties of a given compound. 

In order to generalize the Ginzburg-Landau formula- 
tion to study the nonlinear optical properties in antifer- 
romagnets we first note that the (generalized) Ginzburg- 
Landau functional has to obey the symmetry of the high- 
temperature phase, as the (spontaneous) breaking of this 
symmetry is inherent in the solution which is a minimum 



of the Ginzburg-Landau functional. As an illustration, 
we consider a hypothetical magnetic dipole contribution 
to SHG which is only present in the ordered phase. Let 
us assume the following symmetries for the order param- 
eter, i) it is a c-tensor, i.e. antisymmetric under the time 
reversal operation and ii) it is a (real) axial tensor of first 
rank ( a pseudovector), i.e. Oi. In this case 



F ^ - [x^Jkl{T > Tn) H* E, Ek 
+ x:,kiiT > Tn) H, E; El] Oi 



(6) 



would be a valid expression for the combined free-energy, 
where Xijki{T > T/v) is the susceptibility tensor in the 
paramagnetic phase. The expression in the parenthesis 
indicates that the susceptibility is above the transition 
temperature T^. It should be noted here that the free 
energy is a i-scalar, i.e., a scalar invariant under the time 
reversal operation. This implies that the susceptibility 
tensor in the paramagnetic phase has to be a polar i- 
tensor of fourth rank. Now, comparing Eq. (o) with Eq. 
(g), one would obtain the following relation between the 
SHG susceptibility and the order parameter: 



xgr(T < Tr,) 



Xi]ki{T > Tn) Oi 



(7) 



It is clear from the above expression that from the 
knowledge of the order parameter and the symmetries of 
the susceptibility tensor [ |lO[ in the paramagnetic phase, 
it is possible to obtain all the non-zero components of the 
SHG susceptibility in the ordered phase. Moreover, the 
SHG susceptibility tensor below the transition tempera- 
ture becomes directly proportional to the order parame- 
ter, which will ultimately manifest in the non-reciprocal 
optical effects in the system. In what follows, we examine 
this in detail for the cases of the magnetoelectric com- 
pound Cr203 and of the ferroelectric-antiferromagnetic 
material YMnOs. 

Cr203, in its paramagnetic phase (above T^v ~ 307 K), 
crystallizes in the ccntrosymmetric point group 3m. The 
unit cell contains four Cr"^+ ions, which occupy equiva- 
lent c-positions along the optical axis. This structure has 
a centre of inversion and parity considerations allow only 
axial i-tensors of odd rank and polar i-tensors of even 
rank. Thus, above Tn, SHG electric dipole effects are 
forbidden but magnetic dipole effects are allowed. Below 
T/v, the four spins in the unit cell order along the op- 
tical axis in a non-centrosymmetric AFM structure (the 
spin order being, up, down, up, down) which leads to two 
types of domains transformed into each other by time re- 
versal symmetry. Both, space and time reversal symme- 
try operations are separately broken by the spin ordering 
but the combination of both the symmetries remains a 
symmetry of the crystal. The magnetic point group of 
Cr203 being 3m allows new tensors, i.e., polar c-tensors 
of odd rank and axial c-tensors of even rank. Thus, elec- 
tric dipole effects due to polar c-tensors of odd rank are 
allowed below Tj^. 

In order to construct the nonlinear free energy for 
SHG, one needs to know the order parameter in Cr203. 



This can be achieved following Dzialoshinskii's | [Tl| ] 
derivation for the Landau theory of second order phase 
transitions Q. We assign a mean spin S^ (/3 = 1, 2, 3, 4) 
to each of the four Cr'^^ ions. Out of the Sp one can 
form the following linear combinations m, li, I2 and I3: 



m = Si + S2 + S3 + S4 
li = Si — S2 — S3 + S4 

12 = Si — S2 + S3 — S4 

13 = Si + S2 — S3 — S4 



(8) 



The components of the vectors m and 1^ {a — 1,2,3) 
can be classified according to the irreducible representa- 
tions of the paramagnetic group 3m as given in Table || 
iy,y. The components m^, liz, I2Z, hz transform ac- 
cording to the one-dimensional representations A2g, Aig, 
A-iu, ^2m of the point group 3m. The x and y compo- 
nents of the vectors m and li transform according to the 
two-dimensional representation Eg-, and the x and y com- 
ponents of the vectors I2 and I3 transform according to 
Eu- We make now the following key observation: 

TJie order parameter of a magnet is given by that irre- 
ducible representation of the paramagnetic point group 
which is invariant under the symmetries of the magnetic 
group. 

This statement can be considered as a generalization of 
von Neumann's principle (see pC|) and follows from the 
observation that the thermodynamic expectation value 
of any combination of the constituent spins S^ which is 
not invariant under the magnetic group, i.e. 3to for the 
case of Cr203, can be shown to vanish identically. For 
Cr203 the irreducible representation which is invariant 
under im is Aiu, i.e. hz. 

Thus, the order parameter for Cr203, which is the 
staggered magnetization constructed from I2Z is a c-axial 
scalar. The description given by ^22 corresponds to a 
collinear ordering (up, down, up, down) which is the spin 
ordering in Cr203. Now, we can write down the free 
energy due to the electric dipole contribution as. 



F^-[x^MT>TN)E:EJEk 
+ X*MT > Tn) E, E* E*,] hz 



(9) 



and therefore the SHG susceptibility tensor can be ob- 
tained as, 



Xgf ^(T < Tn) = X^3k{T > Tn) hz 



(10) 



i.e. a c-polar tensor of third rank is obtained from the 
tensorial product of an i-axial tensor of third rank with 
a c-axial scalar (or pseudoscalar). Since the suscepti- 
bility tensor Xijk above T/v is an i-axial tensor of third 
rank, we know all the non-zero components from the sym- 
metry analysis Q. Thus, using Eq. ([lO|), one obtains 
all the non-zero components of the SHG x^ which are, 
Xyyy = -Xyxx = "Xxyx = ~Xxxy Furthermore, the 
symmetry of x{T > TV) in the paramagnetic phase dic- 
tates that x^^'^{T < Tn) is hnearly dependent on the 



order parameter, which is the reason why one observes 
the AFM domains through SHG in CraOs. 

The second example we want to illustrate is YMn03. 
YMn03 is a ferroelectromagnetic material whose crys- 
tal structure above the Curie temperature Tc ~ 913K 
is presumably centrosymmetric and described by the 
point group 6/mmm. The elementary unit cell con- 
tains six Mn'^^ ions. Below Tc, YMn03 orders fer- 
roelectrically and the charge ordering breaks inversion 
symmetry. The YMnOs structure is then described by 
the non-centrosymmetric point group 6mm. The vec- 
tor p('5'-P) = (0, 0,P2'- ') of spontanous polarization is 
directed along the six-fold z-axis. The magnetic prop- 
erties of YMn03 arise from the manganese ions Mn''+ 
in the high spin state S = 2. Below the Neel temper- 
ature, T/v « 74:K, the spins of the six Mn'^+ ions in 
the unit cell are ordered antiferromagnetically in a tri- 
angular structure perpendicular to the polar axis. The 
crystallographic, magnetic and electric properties of the 
hexagonal YMnOs and the other rare-earth manganites 
were studied in the sixties and the related data are avail- 
able in ||l^. New data concerning dielectric, magnetic, 
infrared and Raman studies have been also reported re- 
cently |l|,|l|. 

In the electric dipole approximation, SHG is allowed in 
YMn03 below Tc due to the inversion symmetry breaking 
by the ferroelectric ordering of charges ||l^. Thus, from 
the symmetry of the susceptibility tensor in the paraelec- 
tric phase and that of the order parameter which is the 

(SP) 

spontaneous polarization Pz along the six-fold axis, 
one should be able to write down a free energy like Eq. 
(y). 6/mmm is centrosymmetric, which implies that in 
the paraelectric phase only polar tensors of even rank 

(SP) 

and axial tensors of odd rank are allowed. Pz is an 
i-polar first rank tensor (vector) directed along z. There- 
fore, in order to obtain an i-polar third rank tensor in the 
ferroelectric phase which should describe the SHG elec- 
tric dipole contribution, only contractions of the 6/mmm, 
symmetry polar tensors with odd powers of the order pa- 
rameter are allowed. In the lowest order, the free energy 
can be written as. 



F = -[x^MT>Tc)E*E,Ek 
+X*,kz{T>Tc)E,E*El] Pf^) 



(11) 



such that the SHG susceptibility tensor in the ferroelec- 
tric phase is obtained as, 

xgf ^ {Tc>T>Tn)^ X^,kz {T > Tc) P^/'''^ . (12) 

Here the susceptibility tensor Xijki above the Curie tem- 
perature Tc, i.e. in the paraelectric phase, is an i-polar 
tensor of fourth rank. Now, from the symmetry analy- 
sis pffl of this tensor, one gets the non-zero components 
of x(r > Tc) in the paraelectric phase which when con- 
tracted with the order parameter give rise to all the non- 
zero components of x^^'^ in the ferroelectric phase, i.e.. 



Xzzz and Xxxz{^) = Xyyz(3) ((3) denotes all possible per- 
mutations of the 3 indices). Moreover, it follows from 
the above expression that the SHG susceptibility in the 
ferroelectric phase is a linear function of the ferroelectric 
order parameter. 

As already mentioned earlier, YMnOa belongs to the 
crystal class 6mm in the ferroelectric phase. Below Tn, 
the six magnetic ions in the unit cell order antifcrromag- 
netically and perpendicular to the six-fold axis, i.e. three 
spins are arranged in a triangular structure on planes 
normal to the six-fold axis and separated from each other 
by a distance equal to half the lattice period along the 
hexagonal axis. The magnetic ordering in this material is 
non-coUinear but coplanar and can be determined from 
the exchange interactions among the spins only, similar 
to that in Cr203. The corresponding magnetic group is 
6mm. Following Nedlin [O and Pashkevich et al. [Esl, 
we consider the following linear combinations of spins: 



to local planar displacements of the Mn'^+ ions) and form 
the linear combinations p, q, rj and p analogous to s, 1, 
T and (J in Eq. ([L3|). These linear combinations of the 
i-polar vectors also follow the irreducible representations 
of the spatial group 6mm . From the direct product rep- 
resentations of the spin pseudovectors and that of the 
i-polar vectors (which of course is reducible) , one obtains 
the following combination A = cr^ 77^ -I- cr^Tyj*" which be- 
comes invariant under all the symmetry operations of the 
magnetic group. Thus A should be equivalent to -03 in the 
sense that they belong to the same irreducible one dimen- 
sional representation. Furthermore, from the generating 
matrices of the magnetic group, one can also construct 
an invariant c- axial quantity in lowest order, which in the 
present case is a tensor of rank three, R3 (in 6mm, all the 
c-axial tensors of rank smaller than three vanish). Now 
expanding the free energy in terms of R3 one obtains in 
lowest order, 



n 



CTl 



s = Si 

Si — UJ*S2 
Si — UJ*S2 



- S2 + S3 + S4 + S5 + Se 
S2 + S3 — S4 — S5 — Se 

- L0S3 — S4 + (jJ*S5 + LoSq 

* 
T2 = -Tl 

LuSs + S4 — (JJ*S5 — LoSq 

dl ~ CTi* , 



(13) 



where uj* is the complex conjugate of the phase factor 
io = e*'^''^. The magnetic irreducible representations are 
given by some linear combinations 1/' of the spin compo- 
nents (see Table |ll[) . Since the spin ordering is coplanar, 
■0 may easily be expressed in terms of the components of 
the vectors s, 1, t and a written in the cyclic coordinate 
frame |]l^iH as follows. 



s = e 



(14) 



where s^ ^ s^ :^ is^ and similarly for 1, t and a. Here 
the Ci 's are the unit vectors along the axes of the crystal 
coordinate frame (the z-axis coincides with the six-fold 
axis of the crystal lattice). 

From Table O we learn that "01 , • • • , '!/'4 transform ac- 
cording to the one-dimensional representations Ai, A2, 
Bi and B2 whereas 05 and 06 transform according to the 
two-dimensional irreducible representations Ei and E2 of 
the paramagnetic and ferroelectric group 6mm. The ir- 
reducible representation which remains invariant under 
all symmetry elements of the magnetic group 6mm is 
Bi, therefore -03 is a good candidate to be defined as the 
antiferromagnetic order parameter ||l^ for YMnOs. 

■03 is a complicated combination of the spin compo- 
nents of the six ions in the unit cell. Let us build a more 
intuitive object which belongs to the same Si irreducible 
representation and where not only spin components but 
also vector components are introduced, whose physical 
interpretation may be that of a local field on each Mn'^+ 
ion. Thus, similar to the spins Si ■ • ■ Se, one can intro- 
duce i-polar vectors Vi • • • Ve (which might correspond 



F — ~ [Xijklmn{T > Tn) E* Ej E]~ 
+ XijklmniT > Tn) Ei E* E'l ] Rimn 



(15) 



and thus the SHG susceptibility tensor can be obtained 
as, 



Xijk (T < Tn) — Xijklmn{T > Tjn) Rln 



(16) 



It should be noted that the c-axial tensor R^ has only one 
independent component, i.,e., Ryyy — —Rxxy{^) which 
can be identified as the order parameter discussed above. 
Performing the sum over Z, m and n in Eq. (Hq) explicitly 
we then obtain an alternative expression for the allowed 
matrix elements contributing to SHG in the antiferro- 
magnetic phase as 



xf,f'^(T < T, 



N 



(17) 



[XijkyyyiT > Tn) — XijkxxyiT > Tn) 
'Xijkxyx{T > Tn) — Xijkyxx{T > Tn)] 1p3 



Thus, from the symmetry properties of the sixth rank 
i-axial tensor Xijkimn above Tn ||19|], together with Eq. 
(O), we get all the non-zero components of the SHG 
susceptibility tensor for T < Tn, which are, xtxx' — 
-Xy^ = -Xy"y'^ = ^Xx^y^ (as onc can easily verify, 
see PI). 

Since YMn03 is characterized by two order parame- 
ters, one for the paraelectric-ferroelectric transition and 
the other for the paramagnetic-antiferromagnetic tran- 
sition, it is natural to ask about the coupling between 
them. We observe here that the SHG tensor in the anti- 
ferromagnetic phase, as given by Eq. (|l7|), is directly pro- 
portional to the even-rank i-axial tensor Xijkimn which is 
not allowed in the high-temperature group 6/mmm for 
T > Tc, which is centrosymmetric. Xijkimn need there- 
fore to be proportional to the ferroelectric order param- 
eter (compare Eq. (O) and Eq. dlT 



X^3klr,^n{T > Tn) = XyfeZmn. (T > T,) P]^^) , (18) 



where Xijkimnz is an i-axial tensor of rank seven which is 
allowed in Q/mmm. Comparing Eq. ( |l8[ ) with Eq. (Ill 
we find immediately: 



supported by the Deutsche Forschungsgemeinschaft, the 
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x|r(r <Tn)^ 



\- - — iv/ Xijklninzy-^ -^ -^c) ^l^ni,n ^z — 

[\ijkyyyz\-^ ^ -^c) \ijkxxyz\-^ -^ ^c) 
-Xijkxyxz{T > Tc) - XijkyxxziT > Tc) ] V'S Pz^^'' ■ 

It is therefore possible to derive all the components of the 
SHG susceptibility tensor below T^ from the non-zero 
components of Xijkimno together with both the order pa- 
rameters. It is also clear from the above equation that the 
SHG susceptibility below T^ is directly proportional to 
the bilinear combination of both order parameters which 
in principle, could be verified from experiments. 

The method described in the present manuscript to 
study the nonreciprocal optical effects in the magnetic 
materials through SHG is purely phenomenological and 
based on symmetry considerations. Therefore, it should 
be possible to generalize this phenomenology to other 
hexagonal rare-earth manganites RMnOa where R=Ho, 
Er, Tm, Yb, Lu etc. 

In conclusion, we summarize the main findings of the 
present paper. A phenomenological Ginzburg-Landau 
theory has been developed for second harmonic gener- 
ation in materials which undergo one or more phase 
transitions by expanding the free energy in terms of the 
order parameter/s and the susceptibility tensor in the 
high-temperature phase. We have shown how to ob- 
tain explicitly the SHG susceptibility components as a 
function of certain susceptibility tensors allowed in the 
high-temperature phase and of the order-parameter. We 
have carried through this prescription for the magneto- 
electric compound Gr203 as well as for the ferroelectric- 
antiferromagnetic material YMnOa. We also argue that 
this analysis can be extended to the other hexagonal rare- 
earth manganites. 
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Irreducible Representations 


A2g 


rriz 


Ma 


hz 


Am 


hz 


A-2U 


hz 


Eg 


f m^,my \ 

\ hxjly J 


Eu 


fl2xj2y\ 
\ hx.hy ) 



TABLE I. Classification of the components of the vectors 
m, /, in accordance with the irreducible representations of the 
crystallographic point group 3m,. 





Irreducible Representations 


Ai 


V'i = -(Tr+r2+) 


A2 


^2 = {t- - r^) 

i)'^^s^ 


Bi 


^3 = (^r + fT2+) 


B2 


V'4 = (-erf +cr+) 


El 


(V'5,7 = S-^Ah,U == S-) 


E2 


(V6,7 = ;+,V'6,77 = r) 
(-06,777 = -T2 ,tp6,IV = Tj^) 



TABLE IL Classification of the components of the vectors 
s, I, T and a in accordance with the irreducible representations 
of the crystallographic point group 6mm. 



